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Abstract 



We present results for the renormalized quartic self-coupling \r and the Yukawa coupling 
Ur in a lattice fermion-Higgs model with two SU(2)i doublets, mostly for large values of 
the bare couplings. One-component ('reduced') staggered fermions are used in a numerical 
simulation with the Hybrid Monte Carlo algorithm. The fermion and Higgs masses and the 
renormalized scalar field expectation value are computed on L 3 24 lattices where L ranges 
from 6 to 16. In the scaling region these quantities are found to have a 1/L 2 dependence, 
which is used to determine their values in the infinite volume limit. We then calculate the 
Ur and \r from their tree level definitions in terms of the masses and renormalized scalar 
field expectation value, extrapolated to infinite volume. The scalar field propagators can 
be described for momenta up to the cut-off by one fermion loop renormalized perturbation 
theory and the results for X R and %)r come out to be close to the tree level unitarity bounds. 
There are no signs that are in contradiction with the triviality of the Yukawa and quartic 
self-coupling. 
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1 Introduction 



The experiments in high energy physics are so far consistent with the predictions of the per- 
turbative Standard Model. Two particles still await experimental confirmation, the Higgs 
boson and the top quark. The top quark seems to be on the verge of discovery. Radiative 
corrections and high precision experiments restrict the mass of the top quark to an interval 
of about 120-180 GeV [0, whereas the experiments of the CDF collaboration give already a 
lower bound of 91 GeV with 95% confidence level 0. This indicates that the Yukawa cou- 
pling of the top quark will be relatively small and can still be treated within perturbation 
theory. Similar restrictions on the mass of the Higgs particle are rather weak ||. A lower 
bound on the Higgs boson mass was obtained in ref. ||], where it was shown that it cannot be 
lighter than approximately 50-100 GeV, provided the mass of the top quark lies in the range 
100-200 GeV. The determination of upper bounds on the Higgs mass is a non-perturbative 
problem and lies outside the scope of weak coupling expansions 0. For recent analytical 
studies going beyond ordinary perturbation theory using large N techniques see ref. [|J. 

It is not excluded by the present status of the experiments that there exists also a fourth 
fermion generation with a neutrino heavier than half the Z boson mass. Since the systematics 
of the masses of the known fermions are unclear until now there is no reliable extrapola- 
tion to the masses of a fourth generation. Calculations of the radiative corrections to the 
p parameter indicate that the mass split within the isospin doublets of a possible fourth 
generation must be extremely small. However, the constraints obtained by a comparison of 
experimental data with perturbative calculations may be doubtful if the relevant couplings 
are large. 

Experimentally it is important to know the largest possible values of the Higgs boson and 
heavy fermion masses m# and irtp in the Standard Model. Theoretically this is an issue of 
self-consistency of the model. Large ratios trh/vr and mp/vR (where vr = 246 GeV is the 
electroweak scale) imply a large renormalized scalar self-coupling Xr and Yukawa coupling 
Ur. But it is well known that the Standard Model is suspected to have the property called 
'triviality' and that Xr and i/r can only increase if the regularization scale A decreases. For 
energies larger than A the theory is not consistent anymore and 'new physics' has to take 
over, limn and/or mp get so large that A gets comparable in magnitude, the theory looses 
self-consistency and new physics is noticeable at that mass scale. This scale is somewhat 
fuzzy because the definition of A depends on the method of regularization and because it is 
not clear how large the ratio m^/A can be before new physics takes over. 

One can get an impression of the sensitivity to A by assuming that the one-loop pertur- 
bative /3-functions are valid for momenta up to A and identifying the bare couplings with the 
running couplings evaluated at A. Integration of the renormalization group equations leads 
to a relation between the bare couplings at the regulator scale and renormalized couplings at 
the electroweak scale. The renormalized couplings vanish logarithmically when A increases 
to infinity at fixed bare couplings. Varying the bare couplings at fixed A the renormalized 
couplings go through an allowed region in the coupling constant space. Examples for such 
allowed regions in a fermion-Higgs model with two isospin doublets and gauge couplings 
switched off are shown in fig. [I], where we plotted the ratio m a /vR_ = y/2\ R as a function 
of ur = mpjvR (we shall use ra# and m a interchangeably). The solid and dotted curves 
represent the boundaries of such allowed regions for two different cut-off values. The cut- 
off for the dotted curve is by a factor two larger than for the solid curve. Such plots were 



first presented in ref. [0] were A was chosen at the GUT scale, and discussed further in ref. ||. 

The upper bound curve for \r does not depend very much on y^ and the lower bound 
curve for \r increases rapidly when the value of y^ is raised. The upper (lower) bound 
curve corresponds to infinite (zero) bare quartic coupling and a bare Yukawa coupling rang- 
ing from zero to infinity. Both curves join in a point where yu and Xr have their largest 
values. When increasing A to infinity, the allowed region shrinks to the origin and its bound- 
ary depends only weakly on A. For smaller regularization scales the sensitivity to A increases. 



The above exploration ignores of course the fact that infinite bare couplings are outside 
the region of validity of the one-loop /3-functions. To overcome this we perform numerical 
simulations using the lattice regularization with lattice distance a. The regularization scale 
may be identified with the largest possible momentum, which is n/a for bosons and n/2a 
for fermions, or simply with the inverse lattice distance 1/a. The three values reflect the 
unavoidable arbitrariness in the definition of A. We shall use A = n/2a as the scale where 
the model breaks down and 'new physics' takes over. In order to retain some scaling we shall 
limit the ratios m^/f/A^O.4, which corresponds to More precise definitions of 

upper bounds have been studied in terms of limits to scaling violations 0, [TO, O, [12], [13 . 



It is reasonable first to investigate the question of upper bounds within a simplified 
fermion-Higgs model and to ignore all couplings in the electroweak theory which are small 
and can be treated well within perturbation theory. For an overview about the recent 
progress in this field using the lattice regularization we refer to ref. |L4 



Using the lattice regularization for a quantum field theory that involves fermionic field 
variables one is confronted with the phenomenon of species doubling: On a hypercubic lattice 
each fermion is accompanied by 15 doubler fermions with degenerate mass. This fermion 
doubling gives rise to the following two problems: a) Eight of these species doublers cou- 
ple with an opposite chiral charge |T3| and spoil the chiral couplings to the SU(2) i ®U(l)y 
gauge fields. Several proposals have been developed to overcome this problem HT4[| , but so 
far none of them has been shown to succeed in formulating a chiral gauge theory on the 
lattice. Recently, a thorough investigation of one of the proposals [|Tj| has raised strong 
doubts on its success [|D| although this is not generally accepted [[TJJ. b) A second and 
more serious problem is caused by the large number of mass degenerate fermion species. A 
straightforward, non-chiral transcription of a lattice fermion-Higgs model with one SU(2) 
doublet on the lattice would lead to a theory with 16 mass degenerate doublets, which is not 
observed in nature. It is therefore important to develop methods which allow to reduce this 
number of fermions. At the moment there exist two different approaches, the mirror fermion 

II, which enable us to reduce the 



method [19] and the staggered fermion approach [20, 21 



number of isospin doublets to one. In a numerical simulation with the Hybrid Monte Carlo 
algorithm (HMCA) this number has to be doubled in both cases. 



In the mirror fermion model |19j each fermion field is paired up with a mirror fermion 
field. This 'doubling' of the fermionic degrees of freedom allows for chirally invariant Wil- 
son mass terms to remove species doublers of the original and the mirror fermions from the 
spectrum The physical mirror fermion can also be decoupled from the spectrum (in 

the fermion-Higgs model) by a proper choice of the bare coupling parameters and the model 



then describes one mass degenerate doublet in the scaling region [£J, |25 



Figure 1: The m a -mp plot for the fermion-Higgs model with two mass degenerate doublets. 
The masses are given in units of v R . The various points in this figure represent our non- 
perturbative numerical results for the ratios m a /vR and m F /v R in the infinite volume limit. 
The values for clvr = irv R /2A, given beside the frame indicate the dependence of the cut-off 
A. The symbols correspond to those in the phase diagram, fig. 2. The vertical and horizontal 
dashed lines give the tree level unitarity bounds for the Yukawa and quartic self- coupling. 
The solid and dotted lines were obtained by integrating the one-loop ^-functions. The cut-off 
parameter for the dotted line is twice as large as for the solid line. 



In this paper we will follow a different proposal |2T], |22| which is based on the 'reduced' 
or 'real' staggered fermion formalism [26. pTfl . The fermion-Higgs model we will investigate 
in this paper has already been explored in ref. |2Sj|, where we presented preliminary results. 



The usual euclidean staggered fermions on a four-dimensional hypercubic lattice describe 
four flavors in the scaling region. By using the 'reduced' staggered formalism the number 
of staggered flavors can be reduced to two. These two staggered flavors are coupled to the 
Higgs field, leading to a model with only one doublet in the scaling region. Since the stag- 
gered flavor symmetry group is discrete, the full 0(4) symmetry is broken by the Yukawa 
coupling to a discrete subgroup and it is necessary to add two scalar field counterterms to 



recover the full symmetry in the scaling region pq| . We will not add these counterterms, but 
present several analytical and numerical results which show that the effects of the symmetry 
breaking are small in the parameter region we are interested in. 



The mirror fermion model and our model with reduced staggered fermions are expected 
to reproduce in the scaling region the target model, which is described by the continuum 
action 

S = 3h + 3 f (1.1) 



with 



S 



H 



d 4 x 
d 4 x 



\ |Tr {(yt(^)} + ^ ±Tr {^} + *° |Tr \{f 



3=1 3=1 



[1.2) 
1.3) 



Here 4>(x) is a 4-component scalar field in the 2x2 matrix notation, ipj(x) is an SU(2) 
doublet and denotes the number of SU(2) doublets. The bare coupling parameters m^, 
Ao and yo are respectively the mass parameter, the quartic self-coupling of the scalar field 
and the Yukawa coupling. 



We can regularize this model by introducing a four- dimensional hypercubic lattice with 
lattice spacing a and replacing the derivatives d^{x) by ((j) x +p, — 4>x)/ a - It is convenient to 
rescale the fields 

<j>{x) -> V2^ <p x /a , if)(x) -> ^/a 3/2 (1.4) 
and reparametrize the coupling parameters 



(amof 



1 - 2A 



K 



!J0 



y 



'Ik 



;i.5) 



As usual we shall mostly use lattice units, i. e. a = 1. The actions ( |1.2| ) and ( |1.3| ) then get 
replaced by 



<f>l<f> x + \(4<f> x -iY 



;i.6) 



N D f _ 

EE \Y1 1(^-7^ 

x 7=1 t fi 



;i.7) 



This model was extensively investigated in ref. |2^, ^0|. Because of species doubling it de- 
scribes 16 x No doublets in the continuum limit (naive fermions). 



Figure 2: Phase diagram of the reduced staggered fermion model with No = 2. The solid 
lines indicate the position of the phase transitions between the FM, PM, AM and FI phases. 
At the points marked by the various symbols, we have performed numerical simulations on a 
sequence of different lattices ranging in size from 6 3 24 to 16 3 24. These symbols will be used 
throughout most of the figures of this paper. The dots represent some points at k — where 
we carried out calculations only on a 12 3 24 lattice. 



In this paper we will use the reduced staggered fermion method which reduces this large 
number of SU(2) doublets to No- The reduced staggered fermion formalism will be recalled 
in sect. 2. For the study of the largest possible quartic coupling we consider only A = oo 
which corresponds to radially frozen Higgs fields, i. e. <f>\.<f> x = 1. 



In the last part of this introduction we shall summarize the main results. The phase 



diagram for the model with = 2 obtained in ref. |28[] is shown in fig. |2|. There are four 
different phases in the K-y plane: A broken or ferromagnetic (FM), a symmetric or param- 
agnetic (PM) phase, an antiferromagnetic (AM) phase and a ferrimagnetic (FI) phase. We 
will restrict ourselves in this paper to the broken or ferromagnetic (FM) phase. The AM 
and FI phases have to be regarded most probably as lattice artefacts since their existence 
seems to be related strongly to the hypercubic lattice geometry. An interesting aspect of 
the phase diagram is that the FM phase extends into the negative k region which is not 



accessible in the continuum parametrization ( |1.5|) . We note that yo — > oo for k \ and 
according to eq. ( |1.5|) yo would become imaginary for k < 0. Our numerical results indicate 
that the whole PM-FM phase transition line falls into the same universality class and there 
is no reason to ignore the negative k region. On the other hand our quantitative results for 
the renormalized Yukawa coupling suggest that the region with k < may not be relevant 
since y R saturates when k \ while keeping vr constant and does not increase any more 
when k is lowered beyond k = 0. The various symbols in fig. |2| mark the points in the 
phase diagram where we have carried out the numerical simulation. The points are in three 
different regions which are labeled by the Roman numerals (I), (II) and (III). 

In fig. |I] we have displayed the Higgs mass m CT as a function of the fermion mass mp. 
The masses on the axis are given in units of the electroweak scale vr ~ 246GeV. The plot 
contains only our infinite volume results for the ratios m a /vR = \/2\ R and rrip/vR = y^. 
The symbols in this figure label the position of the various points in the phase diagram and 
correspond to those in fig. |2|. The values of av r listed beside the figure indicate the cut-off 
values of the various points. The arrows mark the tree level unitarity bounds for the Yukawa 
and quartic self-coupling. The graph shows that the numerical results are very close to these 
bounds which implies that the renormalized couplings are not very strong. Also our other 
results are in accordance with the triviality of the Yukawa and quartic self-coupling. It is 
remarkable that the solid curve, which was obtained by integrating the one-loop /3-function 
is not very far from the numerical results. 

The outline of the paper is as follows: In sect. 2 we introduce the model. Sect. 3 contains 
a perturbative one-loop calculation for the Goldstone and Higgs particle propagators and a 
discussion of the effects of the 0(4) symmetry breaking. In sect. 4 we describe the numerical 
methods which we have used for the determination of the Higgs mass, the wave-function 
renormalization constant of the Goldstone propagator and the fermion mass. Sect. 5 deals 
with the extrapolation to infinite volume and discusses the infinite volume results for the 
renormalized Yukawa and quartic self-coupling. A summary of our results is given in sect. 6. 



2 The model 

We start this section by recalling the reduced staggered fermion method which allows a re- 
duction of the number of 16 fermion species by a factor eight. In the second part we will 
explain how the two staggered flavors can be coupled to the Higgs field. 



Let us start from the naive euclidean action 



(2.1) 



XfJ 



where the field ip x is a usual four-component Dirac spinor on the lattice. Because of the 
fermion doubling phenomenon the action (|2.1|) leads to a model which can be represented in 
the continuum by the following action 



S 



K 



Np 

d 4 x ^2 ipj (x)<j)ipj (x) 

3=1 



N F = 16 



(2.2) 



and it describes 16 massless fermions rather than one. Using the usual staggered fermions 
this number can be reduced by a factor four. This is achieved by performing a spin- 



diagonalization transformation on the fields ip x and ip x in eq. (|2.1| ) 

r x = E(7 x ) a/3 xi , C = E ((t^)* , (2.3) 

/3=1 /3=1 

where 7 X = 7i 1 72 2 73 3 74 4 - After this transformation the action ( [2.1|) is a sum of four identical 
terms 

4 

& = -EE hAxZxt+a - xt+^Xx) • (2-4) 

a=l xfi 

The sign factor r\^ x in ( |2.4| ) is given by r/^ = (— l^i+'-'+^-i and the staggered fields Xx 
and x x are four- component complex Grassmann variables. Using only one component and 
dropping the superfix a, it can be shown that the staggered fermion field Xx describes four 
Dirac fermions in the continuum limit. The flavor and spin indices of the four staggered 
fermions are spread out over the lattice and do not appear in an explicit form. To have 
control over the spin-flavor structure it is convenient to introduce the 4x4 matrix fields |20] 

** = §E l X+b Xx +b , % = | E(7* +6 fe +fe , (2.5) 
b b 

where in contrast to eq. ( |2.3|) we sum over the 16 corners of a unit lattice hypercube, = 0, 1 
and the fields Xx are one-component complex Grassmann variables. The row (column) matrix 
index of the *f> x field represents the spin (flavor) label and vice versa for Since the \l/ 
field contains 16 times as many degrees of freedom as the x field, not all components 
are independent. Their Fourier modes ^(p) defined in the restricted momentum interval 
— 7r/2 < < 7r/2, however, are independent [^]. The staggered fermion action may now 
be written in the form 

Sk = ~Y, iTr(^ 7 ^ x+A - * 2+A 7^) • (2-6) 

Xfl 



This form reduces in the classical continuum limit to the action (2.2) with Np = 4. 

This number of staggered flavors can be reduced once more by a factor two by defining 
the Xx fields on the odd sites and the x x on the even sites of the hypercubic lattice, i. e. 

Xx -> 1(1 - e x )xx , Xx -» i(l + e x)Xx (2.7) 
where e x = (— \y x ^+ x ^+ x 3+xi _ ^he insertion of these relations into eq. (2.5) gives 

^x = \Y.^ + % l -^)Xx +b , f, = lE(7 a4 *) t 2( 1 + ^)x*+». (2-8) 

b b 



When inserting the relations ( |2.8[ ) into eq. ( |2.6| ) one can reproduce the action for reduced 
('real' or 'Majorana-like') staggered fermions 



S K = -3 E V^xXxXx+fi ■ (2.9) 



X/J. 



The restriction of the fields x an d X to odd and even sites corresponds to the projections 

-> 75^75) , + 75^75) ■ (2.10) 



This implies the following structure for the matrix fields 



where if) L , if) R , ipL and ipR are 2x2 matrices. The row and column indices of the ipi and i/jr 
fields are respectively the Weyl-spinor and flavor labels, and vice versa for ip L and ip R . 



The model in eq. ( |2.9| ) is invariant under the staggered fermion (SF) symmetry group 
which includes: 

a) Shifts by one lattice distance 

Xx > CpxXx+p j (2-12) 

where = (-1) X »+ 1+ - +X4 . 

b) 90° rotation 

Xx -> S r (R~ 1 x)xr-i x , (2.13) 

where R = R ap is the rotation x p — > x a , x a — > —x p , x T — > x T with r ^ p,a and Sr(x) = 

c) Lattice parity 

Xx -> (-l) Xl+X2+X3 Xix , (2.14) 

where I — I s is the lattice parity transformation for r = 1, 2, 3. 

d) Global U(l) symmetry 

Xx - e^-xx , (2-15) 

where a is a real phase. 

The transformation (|2.12|) can be interpreted as a discrete flavor transformation |27|, [32|. The 



invariance of ( |2.9| ) under the symmetry ( |2.15| ) implies fermion number conservation. Next 
we shall couple the two reduced staggered flavors to the Higgs field such that we recover in 
the scaling region the target model of eq. ([Lip with one SU(2) doublet. We demand that 



the final form of the action is invariant under the SF symmetry group transformations since 
this ensures the staggered flavor interpretation in the scaling region. 

In order to couple the reduced staggered fermion flavors to the Higgs field we first intro- 
duce the 4x4 matrix 

*=(° })=-Ew- (2.16) 



where tp^ \i = 1, . . . , 4 denote the usual 0(4) components of the Higgs fields. When using 
the relations Q2.11| ) and ( 2.16Q one can show that the action 



£ iTr(^ 7/ ^ x+A - tts+^tt*) + yTr(^^C) 



(2.17) 



reduces in the classical continuum limit to the action of the target model in eq. ( |1.3| ) with 



Nd = 1. After inserting the transformations (|2.8|) into eq. (|2.17| ) the final form of the 
fermionic action in terms of the x fields reads 

Sf = XxXx+ii{Vvx + y£xCv.xWp,x) = -5 J2 XxM xyXy , (2.18) 

xp. x,y 

where ^ 

Vpx = T7Y,¥^-b (2-19) 

iD b 



is the average of the scalar field over a lattice hypercube. The hypercubic Yukawa coupling 
arises naturally from the summation over the corners of a lattice hypercube which we intro- 
duced in the definition of the \I/ matrices in eq. (2.8). The fermion matrix M in eq. (|2.18| ) is 



antisymmetric and real. The final form of the lattice action at A = oo is then given by 

S = 2kJ2 <Pia<Pvnx+(l XxXx+^Vfix + ySxCfixVux) 7 ( 2 - 20 ) 



Xfl Xfl 



with Y/fi=i ¥> 2 fix = 1- The action ( |2.20|) is invariant under the SF symmetry group if the Higgs 



field transforms in the following way under: 

a) Shifts by one lattice distance: 

W C 1 ~~ 2 <W)W,*+,3 , (2.21) 

b) 90° rotations: 



where n = (|, |, \, \ 



ffix Rfiu L Pu,R-^{x+n)-n , (2.22) 

■ 

c) and lattice parity: 

W -> (25^4 - tyVnJx ■ (2.23) 
The additional shifts by the vector n in eq. ( |2.22j ) are needed because the relation between 



the \l/ and x fields in eq. ( |2.8| ) is not manifestly rotationally covariant. The invariance under 
rotation would have been more transparent if we would have associated the \1/ and p fields 
from the beginning with the dual lattice which is shifted by the vector n with respect to the 
lattice for the x fields. The scalar field does not transform under the U(l) symmetry ( 2.15| ). 



The action is not invariant, however, under the full 0(4) flavor group, but one expects 
to be able to recover this invariance in the scaling region. In the scaling region operators 
with dimension larger than four become irrelevant. There are, however, two operators with 
dimension four which respect the discrete symmetries (|2.21|) -( p.23 ), but break 0(4): 



= E <P%, = \ $>^ A - Vfix) 2 - (2.24) 

Xfl ^ Xfl 

We will show in the next section that these terms are indeed generated by the quantum 
fluctuations. In order to recover the full 0(4) symmetry one has to add these operators as 
counterterms to the action (|2.20|) 

S -> S + e O {1) + 5 O {2) , (2.25) 

and tune the coefficients Eo and Sq as a function of the bare parameters k and y such that 
the 0(4) invariance gets restored in the scaling region. In this paper we will not add these 
counterterms, but present several analytic and numerical results which show that the effect 
of the symmetry breaking is small in the parameter region of interest. 



3 One-loop fermion effects on the scalar propagator 

In the pure 0(4) model the numerically measured scalar field propagator has a momentum 
dependence that is nearly of the free field form because the renormalized self-couplings are 
small. The Yukawa interaction with the fermions affects the scalar propagator in three ways: 
the masses and wave-functions are renormalized, the additional self-energy gives rise to a 



more complicated momentum dependence than the almost pure pole found in the 0(4) model, 
and with our staggered fermion method the fermions induce 0(4) symmetry breaking effects. 
The last two effects can be studied in one fermion loop renormalized perturbation theory. 
This is useful even at large bare couplings because the renormalized couplings turn out to 
be relatively small. By taking into account the one-loop effects we can then more reliably 
extract the scalar masses and wave-function renormalization constants from the numerical 
data and estimate the effect of 0(4) symmetry breaking on the renormalized couplings. 

If the fermion effects are neglected the scalar sector of the model has the following 
approximate effective action, 



1» a , m R i \2 



(3.1) 



where d^tp x = ip x +p,~ l Px an d the subscript R denotes renormalized quantities. In the broken 
phase we can decompose the scalar field in a Higgs mode, a, and three Goldstone modes, 
7r a , a = 1, 2, 3, according to 

PR? = v R el + <p' li = (v R + a R )el + 7r£e£ . (3.2) 

The {e^} form an orthogonal set of 0(4) unit vectors and v R is the scalar field expectation 
value. Before taking fermion effects into account, the choice of e^'s is arbitrary and after 
substituting the decomposition (|3.2|) into the action (|3.1j) one finds the usual tree level 
relations for the masses of the Higgs and Goldstone modes, 

2 

ml = 2X R v 2 R , ml = 0, v 2 R = ^, (3.3) 

*R 

as well as the three point interactions \rVr[<j r + o-r7t r 7t r ]. The scalar propagator 

G M = (^T,p»M^-y))) (3.4) 

\ v x,y I 

is given by 

= l^(k 2 +ml)+t^(k 2 + ml), k 2 = 2^(1 - cos (3.5) 

A* 

where and t^ v are the longitudinal and transverse projectors onto the a and tt subspaces. 

To find the effect of the staggered fermions on the scalar propagator we compute the 
one-loop diagrams shown in fig. ||. The Feynman rules for the fermions can be derived from 
the action 

Sf = -5 £ XxXx+p,{Vnx + mpele^x) - \vr £ XxXx+pSp'^Ci^ > ( 3 - 7 ) 

Xfl Xfi 

which results after inserting ( |3.2| ) into eq. ( 2.18 ). The fermion mass in the first term is given 
by the usual tree level relation 

m F = VrVr . (3.8) 
The bar on ip' indicates hypercubic averaging, as in ( 2.19|) . 



Figure 3: Feynman diagrams for the one fermion loop contributions to the vacuum expec- 
tation value (a) and to the scalar field self-energy (b and c) in the FM phase. 



We shall use in the following the staggered fermion formalism developed in refs. |27], [32 . 
The fermion propagator follows from the first term in eq. 



SAB(p)$(p-q) = E e 



-i(p+-KA)x 



(XxXy)e 



i(i+ 7r B)y 



x.y 



Sab{p) 
5(p - q) 



S^H^b sinp^ + mpe^^T^AB cosp M ] 
E^sin 2 ^ + (e^) 2 m 2 F cos 2 ^] 

V5 PA mod(27r) . 



(3.9) 
(3.10) 



The lattice momentum of the fermion is in the restricted interval p^ G (— 7r/2, 7r/2], 7Ta, 
A = 1, . . . , 16, are the momentum four- vectors with components equal to or tt and V is the 
lattice volume. The 16-dimensional gamma and flavor matrices r M and H M form a Clifford 
algebra with {T^, T^} = 25^, {H^E^} = 25^, [T^E, 



0. 



The second term in (3/7) contains the interaction with the scalar fields averaged over a 
hypercube. The corresponding vertex function is given by 



^^ab{p, -q, k) 
h(k) 



-yRh{k)e ikti/2 cos(g 



2 k^je 



ITT B 



l 8(p - q + k + TT A + 7lB + 7T £ + 7TA, 



16 



E 1 



-ikb 



-ikn 



n 



cos • 



(3.11) 



where n = (5, |, 5, 5) is the vector introduced earlier in Q2.22j ). The quantities p + it a and 
q + ttb are the outgoing and incoming wave vectors of the fermion field, with p and q in 
the restricted interval mentioned above, and k is the outgoing momentum of the scalar field 
ip'^. The ti e and are defined such that e x = exp(iii £ x) and = exp^m^x). The factor 
h(k) is due to the hypercubic fermion-scalar coupling. In the classical continuum limit the 
momenta in lattice units approach zero and 



r M (p, -q, k) -> -y K S 5 r 5 E: At 5(p - q + k). 



(3.12) 



Note that yn is a bare Yukawa coupling as it does not contain at this stage the effects of the 
fermion interactions. 



Let us first consider the one loop fermion effect on the vacuum expectation value of the 
scalar field (cf. fig. 0a): 
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ikx 



v R ei5(k) - \N D £ 
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(3.13) 



+ E f % + h) A W 

ml) J 

\y R m F N D el E 



cos 2 Pj, 



E p [sin 2 p p + m 2 F (ef) 2 cos 2 p p ] 



(3.14) 



where iV^ is the number of fermion doublets (which is two in our numerical work) and E P is a 
normalized sum over the lattice momenta. The direction of spontaneous symmetry breaking 
is compatible with the fermion loop correction if the latter has no transverse component, 
i-e. Ei/ = °- This is the case for e t = (±1,0,0,0), . .. , (0,0,0,±1), (±1,±1,0,0)/V2, 
. . . , (±1, ±1, ±1, 0)/y/3, (±1, ±1, ±1, ±l)/2. By studying the one fermion loop effective 
potential 

_J_ Tr 1ti E M Hr M sinp M + yRtpR^EnEsTs cos^] 

2 



Veff(<p) 



16 



£ M [sin Pfl + y 2 R <p 2 R „ cos 2 p M ] 



(3.15) 



for small y R we found in ref. [28| that only = (±1, ±1, ±1, ±l)/2 are local minima which 
correspond to a ground state; the others are saddle points or local maxima. We checked this 
conclusion by a numerical study of the ((p^) probability distribution. In the following we 
shall understand ei to be one of the 16 unit vectors 



1 



±1,±1,±1,±1). 



(3.16) 



Note that the hypercubic coupling prohibits a fermion induced staggered magnetization 
{ipptx) e x . One way to see this is that the factor h(k) in T^(p, —p,k) in ( |3.13| ) vanishes 



when kfj, — it. 



We now turn to the fermion contribution to the scalar self-energy corresponding to the 
diagrams b and c shown in fig. [3]. The tadpole contribution of fig. ^b is given by the following 
expression: 



Eg (*) 
D{p) 



cos pi 
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E 



sin 2 p p + -m% cos 2 p p 



(3.17) 
(3.18) 



which is independent of k. Diagram c in fig. H leads to the contribution 



\y 2 R N D \h(k)\ 2 e l 2^- K) 



cos(p + lk) M cosQo + \k) l 



— Ti{E^T 5 E fl S(p)E 5 T 5 E fl S(p + k)} 



+ \m\ cosp • cos(p + A;)] + Z^m^fcos^ cos(p + k) v + n v]\ . (3.19) 
For k = this reduces to 



Eg(0) = W R K 



D 
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-<W£t5 + ^ 2 "4E D 2 

p U p u 



(3.20) 



where D was denned in (|3.18|) and we use the notation 



c M = cosp^, Sf, = smpp, s 



E 



srn 



(3.21) 



The total fermion contribution to the scalar self energy equals T,^ u (k) = E^(fe) + Ejfj(A;). 
In an 0(4) symmetric model the transverse parts of the zero momentum Ejjj and E^fj in 
( |3.17| ) and ( |3.20|) would cancel such that m^ remains zero. Using 5^ = + we see such 
a cancellation in the first term of ( |3.20| ). However, despite its factor l^, the second term 
also contains a transverse part, as follows from 



c 2 c 2 



c 2 c 2 



l \xv 2-^i JJ2 + V4L 



D 2 



(3.22) 



So we find 



E^(0) 



ylN D U2^ + 2m% C M ) + 2s R v 2 R 



1„,2 
2 



(3.23) 
(3.24) 



leading to a pion mass oc Er. The values of f e range from 0.0054 to 0.0043 for rriF ranging 
from to 0.5. 



In ref. |28| we followed a slightly different strategy and computed the coefficient er of 
the term Er <p R/M m the effective potential fl3.15| ). There we found the same coefficient Er 
as computed here from the two point function. 



The continuum limit of E Atl/ (/c) can be calculated by separating the integration region 

v -> /-!/2 rf4 P/7r 4 ) info 
which leads to the form 
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(3.25) 



Here c_2, Co and r are numerical coefficients which we have not calculated explicitly. The 
term SRS^k 2 corresponds to the second 0(4) symmetry breaking term of dimension four in 



eq. ( |2.24j ). The coefficient Sr can be expressed as 
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[3.26) 



with a = 1 — (The expression ( |3.26|) differs from that in ref. ]28| which was incom- 

plete). The values of f$ range from 0.0026 to 0.0017 for rrip ranging from to 0.5. 



The quadratic and logarithmic divergencies of the continuum limit correspond to the c_2 
and Yn.m 2 F terms respectively. Note in particular that Sr and Er are finite (the continuum 
limit of Er and Sr corresponds to letting nip — > in (|3.24j) and ( |3.26| )) . This is because the 
integrand near the origin in momentum space is covariant, as in the continuum. At the two 
loop level Er and Sr probably acquire logarithmic divergencies with corresponding renormal- 
ization group /3-functions. Presumably this means that Er and Sr will be further suppressed 
by triviality. For example, in pure O(N) models /3-functions for asymmetric couplings have 
been studied previously [j33| , and from these results we found it amusing to discover that 



the ratio Er/Xr in the asymmetric 0(4) model vanishes logarithmically slowly in the scaling 
limit. 



We now perform mass and wave function renormalization. A change in the parameter 
m R m flOD allows for a renormalization of c_2 like terms and a change of scale of the scalar 
field allows for a renormalization of the r-like terms. We may then write for the renormalized 
propagator 

Gll v {k) = 2{\' R + e R )v R l^ + 2E R v 2 R t llu + k 2 S, u + SrS^I + T^\k) , (3.27) 

~ d 



SjT^) = X, v (k) - E M „(0) - k] 



dk 2 ^ 



(3.28) 



k=0 



Here we defined X' R by its appearance in (|3.27|) . This definition is based on the effective tree 
level action 



S, 



eff 



(3.29) 

which leads to fl3~27|) with E( su6 ) -> 0. 



To determine the scalar propagator from (|3.27|) we use the a, it basis, 
Gf(k) = ey u G R ^(k), m b) (k) = ey^\k) 



G a R(k) = IfiuG^lk) 



GnR^k) — ^tfii/G^v — ^ 

a=l L 



2(\'R + E R )v 2 R + (l + 6 -f)k 2 

)k 2 + s aa 



1 -1 



2£rv 2 r 



1 + 



6, 



R) 

-1 



(3.30) 
(3.31) 

(3.32) 



In the cr, 7r basis the transverse and longitudinal projectors are of course diagonal, but 
the Sr term has an off-diagonal piece. To establish ( [3.3 1| ) and ( |3.32|) for the scalar prop- 



agators it is convenient to use the basis vectors e 1 



and e 



-1,-1,1) 



-1,1,-1,1 



\{— 1,-1,1,1) for the modes orthogonal to e 4 = 7,(1,1,1,1). With this choice 
the a, 7r diagonal part of the Sr term takes an 0(4) invariant form, SRS a/3 k 2 /4, and gives 
an additional wave-function renormalization by a factor 1 + Sr/4. The off-diagonal terms 
<5r X^( e ^ e fi — ^ a/3 /4)^ lead to corrections of order S 2 R in G a ^ which may be neglected. 



From eq. ( 3.31 



) and Q3.32 ) we read off the following one fermion loop estimate for the 



Higgs and Goldstone masses, 



, _ 2(\' R + e R )v R , _ 2e R v R 



< - 1+^4 • < - imi ■ (3 ' 33) 



In ref. |28| we found that this result for the Goldstone mass was in reasonable agreement 
with the numerical data, which will be supported by our results in sect. 4. This motivates 
us to take also the above result for the Higgs mass seriously and use it for an improved 
definition of the quartic coupling which is corrected for the 0(4) symmetry breaking, 



We shall compare this with the usual definition 



9 

mi 



A * = 24 (3 ' 35 > 

in our numerical work to be discussed below. After checking that the corrections are small, 
we may neglect them. 

To this one loop order the coupling y R may be interpreted to include also the fermion 
feedback on the scalars. In the spirit of renormalized perturbation theory, ( 3.31 ) and ( p. 32 ) 



are expected to be valid whenever the renormalized couplings are sufficiently small, even 
when the bare couplings are large. 



It can be seen from (|3.25|) and (|3.28|) that becomes important for fc 2 ^4m^. We 



shall use expressions ( p.31| ) and ( |3.32| ) in analyzing our scalar propagator data, using for 
S(swfe) its full lattice form and k 2 — > k 2 as in ( |3.5|) . 

4 Finite volume results for the fermion and scalar prop- 
agators 

In this section we turn to the numerical simulations. After reporting on some technical 
details, we discuss our methods for analyzing the fermion and scalar propagators and we 
present the finite volume results. 

We have measured the propagators in three different regions of the bare (k, y) parameter 
space, k = 0.29-0.31, y = 0.7 (I), k = 0, y = 3.6-4.8 (II) and k = -0.25, y = 5.6-5.8 
(III). These regions with the points at which we have performed the simulations have been 
indicated in the phase diagram of fig. 0. The region (I) was included to get a crude idea 
of the yn dependence of Xr for intermediate values of the Yukawa coupling. The regions 
(II) and (III) were chosen to study the model at large values of the bare Yukawa coupling. 
Region (II) would correspond to yo = oo in the continuum parametrization of the action 
given in eq. ( |1.3|) . 

We have restricted our calculations to the FM phase, in particular, we did not increase 
y and lower k further beyond region (III), since from comparing the results in region (II) 
and (III) we observed that the renormalized couplings \r and remain almost constant 
when lowering k beyond k = 0, while keeping civr roughly fixed. We expect this trend 
to continue also beyond the multicritical point A. The bosonic particle spectrum in the FI 
phase, however, differs from that in the FM phase. It contains in addition to the usual 



Figure 4: The quantity sinp 4 / S^{p^) as a function o/sin 2 p 4 for k = and several values 
of y. The lattice size is 12 3 x 24. The straight lines were obtained by fitting the Ansatz ( \j-4 ) 
to the numerical results. 



Higgs and Goldstone particles so-called 'staggered' Higgs and Goldstone particles which are 
associated with the antiferromagnetic ordering in that phase [M. 

We carried out the simulation on volumes of size V = L 3 x T where the spatial extent 
ranged from L = 6 to 16 and the time extent was always kept fixed at T = 24. The 
trajectory length r in the HMCA was set equal to one and we have tuned the step size St 
such that the acceptance rate fell into a range between 70 and 80%. In general St had to be 
reduced slightly when the value of k was lowered and the value of y was raised. We used for 
the fermions periodic boundary conditions in spatial directions and antiperiodic boundary 
conditions in the time direction. The scalar fields had periodic boundary conditions in all 
directions. Depending on the autocorrelation times for the scalar field expectation value we 
have generated 5,000 to 20,000 trajectories, which resulted in reasonably small statistical 
errors for the various observables. We could afford large lattices and high statistics because 
the staggered fermion matrix in eq. ( 2.18| ) is relatively small and the conjugate gradient 



inversions were found to converge excellently. 
4.1 Fermion propagator 

The fermion propagator in the momentum space representation is defined by the expression 

Sab(p) = (^J^e^ + ^M-h-^ + ^ , (4.1) 

where M xy is the fermion matrix defined in fl2.18j ). Assuming that loop effects are small, we 
can use the free fermion result in eq. (pT9|) to parametrize Sab- This suggests to measure 



the projections 

Sj = i±Tr{T,S} , = ^Tr{r 5 S 5 S,5} . (4.2) 

On a finite volume, the scalar field will tunnel from one minimum of the effective potential to 
another, where the average field value J2x ^^x/V is proportional to (±1, ±1, ±1, ±1) in these 
minima. The propagator component ST is the same in all minima, but is proportional 
to Yjx (ffix/V and would vanish in a finite volume. 

In practise we have measured S^(p) by averaging over rotated scalar field configurations. 
These global 0(4) rotations were chosen such that J^x^ux/V oc 5^ on all configurations 
that went into the averaging. The fermion momenta were chosen along the 4-direction: 

P4=(^--j-, n = - I + l I ..., I (4.3) 

and pj = 0. Typical results are shown in fig. £|, where we have plotted smp^/S^ (p 4 ) for sev- 
eral values of y and k = as a function of sin 2 p4. The linear behavior which is characteristic 
for weakly interacting fermions is seen to hold remarkably well over the full range of momenta. 

In order to extract the fermion mass from these data we adopt the free fermion parametriza- 
tion, 

sinp 4 (1 - m F ) sin 2 p 4 + mp 

Wijhj ~ Z~f ' {AA) 

and fitted the parameters mp and Zp. This formula follows from (|3.9|) with the modifica- 
tion that the vector in the mass part, which corresponds to the direction of symmetry 
breaking, now is the unit vector in the 4-direction. 

The results of the fits are collected in table [3]. The fermion mass found from Sf with 
the appropriate free fermion fit gave the same results within 10%. This small discrepancy 
decreased on larger lattices, presumably because (<£> M ) gets frozen in one of the minima. 



4.2 Scalar propagator 

We have measured the following momentum space Green functions, 

1 3 

G a (p) = G M (p), G n (p) = - G mm (p) , 

d m=l 



(4.5) 



with 



GM 
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xy 



jipix-y) 



(4.6) 



which is identical to (|3.4j ) for non-zero momenta p. The scalar fields in G liu (p) have been 
rotated such that J2x <Px/V oc <5 M 4, as was done for the computation of the fermion propa- 



gator. For the definition of G a ^ in ( |3.30|) - (|3.32| ) this is just a change of the basis such that 



&u * ^V4j G fj Ifiv^fii/ * ^44) »-» 7r 3"/Lti/ v - J <fj,ii ' 3 L-im=\ ^ rami 

that our fitting formulas ( |3.31| ) and ( |3.32| ) should be unchanged, i. e 
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3 ^fiv^* iiv 



\YL=\Gmm, as in (fOP . Note, however, 



we still use in S 



(sub) 



the 



e A of eq. (|3.16| ) (The dynamical fermions in the Monte Carlo process 'experience' of course 



the preferred local minima ( |3.16 ) and not e 
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Table 1: Results for the fermionic observables mp, Zp, and y in the regions (I), (II) 
and (III) of the phase diagram (c.f. fig. [|j for several values of L and T = 24. 



In fig. we show a representative example of the measured inverse propagators G (Tt7r (p)~ 1 , 
as a function of the bosonic lattice version of the momentum squared, p 2 = 2 1^(1 — cosp^). 
To reduce statistical errors we have computed the propagator in the small momentum range 
for all possible lattice momenta, and averaged over results with the same value of p 2 . 



For sufficiently small momenta, p 2 <C m 2 F , one expects a linear p 2 dependence, 

G~\ = (p 2 + ml^/Z^ and a naive linear fit was commonly applied to extract and 
Z a ,n However, the data in fig. |^ show a significant curvature, which may give 



a bias in an estimate of m CT)7r and Z a ^, obtained from a linear fit. This bias can be large 
on small lattices, where even the smallest momenta are too large to neglect the curvature. 



Therefore we followed the same strategy as in ref. |3(J and included the subtracted self- 
energy E( su6 )(p) (cf. eqs. ( |3.31| ) and (|3.32| )) to parametrize the non-linear p 2 dependence. 
Furthermore it was found that the bending over of G~\ for very large momenta, p 2 ^4 (see 
fig. HD, could be reproduced by including the pole contribution coming from the 'staggered' 
scalar modes. Apart from the usual Higgs and three Goldstone bosons the spectrum in the 
FM phase contains also four 'staggered' scalar modes which are associated with the antifer- 
romagnetic ordering in the FI phase. These modes produce a pole in the diagonally opposite 
corner of the Brillouin zone, at 16 — p 2 = (m|*) 2 . The mass m|f of these particles is usually 
large, but it becomes small when approaching point A in the phase diagram. See ref. f3(| 
for a detailed discussion of the staggered scalar modes. 



These considerations lead us to use the following parametrization for the inverse scalar 
propagators in the FM phase, 

G^{p) « Z^/ [f + m 2 ^ + ^£ (su6) (p)^ + Z|7(16 - p 2 + {mtf) . (4.7) 

This requires us to fit five parameters: the mass and wave-function renormalization factors for 
the scalar particles (a or n) and for the staggered bosons and the coefficient of the subtracted 
self-energy, y. This last parameter measures the strength of the renormalized three point 
Yukawa interaction and should be close to yn obtained from the tree level relation with the 
fermion mass. Notice that Yi( su b)/y 2 R contains no free parameters, since we use the measured 
value of the fermion mass in the loop sum. The subtracted self-energy S( snil )(p) is computed 
numerically for the lattice sizes used in the simulations. The coefficient of the term linear in 
p 2 is computed with the finite difference approximation [<9£(fc) / 'dk 2 ] k=0 « (£(/i 4 ) — S(0))//?-4, 
where is the smallest non-zero lattice momentum in the time direction. 



If the renormalized Yukawa coupling is weak, one expects that the one-loop expression 
( |4.7| ) gives a good description of the non-linear p 2 dependence of the inverse propagator. 
Conversely, an accurate description of the data over a wide momentum range with the one- 
loop form (4/7) indicates that the Yukawa coupling must be weak. As will be seen shortly, 
this is the case in our model and including the self-energy term leads to stable values for the 
mass and wave-function renormalization factor even on small lattices. 



We applied this method of analysis to the scalar propagators in all three regions of the 
phase diagram. In region (I) we found larger autocorrelations which urged us to increase the 
statistics to 20,000 trajectories. Region (III) was not very different from region (II). Since 
the AM phase is closer here, we found a stronger effect of the staggered scalar modes on the 
scalar propagator for large momenta, but this could still very well be fitted with the Ansatz 



Figure 5: The inverse scalar propagators G~ 1 (p) (fig. a) and G~ 1 {p) (fig. b) as a function 
of p 2 at the point (n,y) = (0,3.8). The lattice size is 12 3 24. The circles were obtained by 
fitting the Ansatz U-V t° the Monte Carlo data which are represented by the crosses. 



Figure 6: The inverse scalar propagator G n (p) is plotted for the full momentum range as 
a function ofp 2 . The computation was performed at the point (K,y) = (—0.25,5.6) (region 
(III)) on a lattice of size 12 3 24. The points in the lower figure were obtained by fitting the 
Ansatz (\4-Tj ) to the Monte Carlo data which are displayed in the upper figure. 
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12 


0.235(6) 


0.13(4) 


0.65(1) 


1.28(7) 


1.33(2) 


2.75(9) 




0.31 




6 


0.346(7) 


0.13(3) 


0.84(3) 


1.04(4) 


1.20(4) 


2.4(1) 








8 


0.313(8) 


0.10(4) 


0.80(3) 


1.21(6) 


1.23(4) 


2.6(1) 








10 


0.30(1) 


0.14(4) 


0.84(2) 


1.34(9) 


1.28(2) 


2.8(1) 








12 


0.290(7) 


0.14(4) 


0.81(2) 


1.38(7) 


1.27(2) 


2.79(9) 


/tt\ 

(II) 


0.00 


3.6 


12 


0.086(3) 


0.15(2) 


0.28(2) 


0.26(1) 


0.289(3) 


3.2(3) 






3.8 


6 


0.187(4) 


0.15(1) 


0.74(4) 


0.34(1) 


0.56(1) 


4.0(2) 








8 


0.157(3) 


0.16(2) 


0.64(1) 


0.35(1) 


0.500(4) 


A A f A \ 

4.1(1) 








10 


0.143(2) 


0.11(2) 


0.55(2) 


0.35(1) 


r\ A t~i A f v \ 

0.464(5) 


3.9(1) 








12 


0.134(3) 


0.11(3) 


0.50(1) 


0.35(1) 


r\ A A -t / C\\ 

0.441(3) 


3.69(8) 








16 


0.127(3) 


0.17(3) 


0.50(2) 


0.36(1) 


0.410(7) 


3.9(2) 






4.0 


6 


0.203(4) 


0.18(2) 


0.81(3) 


0.40(1) 


0.58(1) 


4.0(2) 








8 


0.185(2) 


0.15(1) 


0.72(1) 


0.412(6) 


0.555(6) 


o r\r\ { pt\ 

3.90(7) 








10 


0.181(2) 


0.12(2) 


0.68(2) 


f\ A A A ( c\ \ 

0.414(9) 


0.537(7) 


3.8(1) 








12 


0.180(3) 


0.11(4) 


0.65(2) 


0.41(1) 


0.528(9) 


3.6(1) 






A 9 


p. 






n 1 A (o\ 
u. L^yz j 


n 7%( a ^ 
u. i oy^) 


n aa (i \ 


n x.%(9\ 
v.ooyz j 


\ ) 








12 


0.214(4) 


0.18(2) 


0.87(2) 


0.47(2) 


0.64(1) 


4.1(1) 






4.8 


12 


0.293(4) 


0.29(2) 


1.34(3) 


0.54(2) 


0.90(2) 


4.6(1) 


(HI) 


-0.25 


5.6 


6 


0.174(2) 


0.17(1) 


0.76(2) 


0.180(4) 


0.293(4) 


4.4(1) 








8 


0.151(2) 


0.11(2) 


0.65(1) 


0.175(4) 


0.268(3) 


4.34(9) 








10 


0.135(2) 


0.17(2) 


0.59(1) 


0.192(4) 


0.251(2) 


4.4(1) 








12 


0.129(1) 


0.13(1) 


0.53(1) 


0.185(2) 


0.237(1) 


4.13(6) 






5.8 


6 


0.182(2) 


0.18(1) 


0.84(3) 


0.209(4) 


0.321(9) 


4.6(2) 








8 


0.167(2) 


0.16(2) 


0.74(2) 


0.207(4) 


0.302(5) 


4.4(1) 








10 


0.163(2) 


0.11(2) 


0.69(2) 


0.205(4) 


0.288(6) 


4.2(2) 








12 


0.157(2) 


0.16(1) 


0.63(1) 


0.213(3) 


0.271(2) 


4.00(7) 



Table 2: Results for the bosonic observables vr, m a>1T , Z C<K and the ratio m a /vR in the 
regions (I), (II) and (III) of the phase diagram at several values of L and T = 24. 



We have fitted the inverse Goldstone propagators over the momentum range < p 2 < 1.3 
and the a particle propagator over the interval 2(1 — cos(27r/T)) < p 2 < 1.7. In these in- 
tervals the effect of the staggered pole in eq. (|4.7|) is negligible and could be dropped. In 
fig. |5| we show an example for a such fit at a point which falls into the region (II). The 
result of the fit is shown by the open circles and it is seen that the agreement with the 
measured data (represented by the crosses) is excellent. In fig. ^|b we show the result of 
a fit over the full momentum range at a point which belongs to the region (III), this time 
including the staggered scalar pole. Again it is seen that the quality of the fit is very good, 
even for large momenta. The resulting values for m 7rj0 -, Z nj(T and y are consistent with those 
which were obtained from the fit over the restricted momentum interval. The fitted values 
for m 2 ^ and Z a ^ are collected in table 0. The fitted values for y can be found in table |I|. 
They are seen to be roughly equal to the y^ values, which gives further support to the consis- 
tency of our analysis and indicates that the renormalized Yukawa coupling is not very strong. 



By using the one fermion loop fit method we achieved a considerable improvement of the 
scalar propagator analysis. In ref. [28] we have used linear fits, but this was possible only 
on the 12 3 24 lattice since on smaller lattices the systematic errors due to the curvature at 
small momenta were unacceptably large. This problem prohibited a study of the finite size 
effects and the extrapolation to infinite volume. The linear fits on the 12 3 24 lattice produced 
slightly larger values of m a and smaller values of Z w . We emphasize that with the fitting 
method used in this paper the systematic error for m a could be reduced by almost an order 
of magnitude. 



The renormalized couplings Xr and y^, listed in the tables [I] and were obtained from 
the tree level relations ( |3.3|) and ( |3.8| ). For the computation of these couplings we have to 
determine the renormalized scalar field expectation value 

vr = "4f= > ( 4 - 8 ) 

V ^7T 

where v is the unrenormalized scalar field expectation value. As we mentioned in sect. 4.1, 
on a finite lattice the system is observed to tunnel from one minimum to another and after 
averaging over many configurations one would find zero for the scalar field expectation value. 
To compensate for this drift it has proven to be useful to rotate the magnetization vector 



for each configuration first into a certain direction, before taking the assemble average [51 
On a finite lattice v can then be defined by the relation 

^dkU' (4-9) 

* X 

As in the pure 0(4) model we use Z„ for the renormalization of the scalar field expectation 
value in eq. (|4.8|), rather than Z a , because the Higgs particle is unstable. 



One should be aware that the fitted value of Z n contains the factor (I+Sr/4) 1 , cf. (pj] 



which implies that v' R = v R {l + 5r/A) x I 2 should be used for the calculation of corrected 
couplings A'^j and y' R from eqs. (|3.8|) and ( |3.34|) , 



V' R = VR(l + S -f)" 2 (4.10) 



Figure 7: The Goldstone mass as a function of y for k = 0. The lattice size is 12 3 24. The 
measured values ofm n are represented by the squares and the one-loop result for obtained 
from eq. ( 3.33j) by the diamonds. 



The ratios R y = (y R — y' R )/y' R and R\ = (\/2\ R — J2X' R ) / \J2\' R give a measure for the 0(4) 
symmetry breaking effects. These ratios were computed numerically from eq. ( |4.10|) and 
( |4.11|) , where the quantities e R and 5 R were determined by means of eqs. ( |3.24|) and ( |3.26| ), 



with the measured values of the fermion mass and y R supplied. For the data listed in tables 
[I] and we find 

\Ry\ < 5% , \R X \ < 6% , (4.12) 

in a parameter region with mp < 0.5 and m a < 0.7. In general the numerical values of \R y \ 
and \R\\ decrease with increasing lattice volume (see also fig. pT0|) . 



The accurate description of the momentum dependence of the inverse scalar propagator 
over a wide momentum range, using one fermion loop renormalized perturbation theory 
indicates that y R is in the perturbative range and it suggests that also our computation of 
the 0(4) symmetry breaking effects should be reliable. This was found to be the case in 
ref. ||28|| , where we compared the measured values of the Goldstone mass with the one-loop 
result ml oc e R of eq. ( |3.33| ). In fig. |7| we repeat this comparison using our high statistics 
data. Again it is seen that the one-loop result gives a reasonable description of the data. 
The increase of for y values approaching the phase transition around y ~ 3.6 is due 
to finite size rounding which is not contained in the one-loop formula. Motivated by this 
result we shall assume that also the one-loop corrections appearing in ( (4.10| ) and ( f4.11p give 
a reasonable account of the remaining symmetry breaking effects. 



5 Infinite volume results for the renormalized couplings 



From the finite volume results for the fermion and Higgs mass we want to obtain estimates 
of the renormalized couplings y R and \r in the infinite volume. Therefore we need to 
extrapolate the results of the previous section to the infinite volume. In subsect. 5.1 we shall 
discuss how these extrapolations have been done. The resulting infinite volume estimates 
for the renormalized couplings are discussed in subsect. 5.2. 



5.1 1/L 2 extrapolation to infinite volume 

If the spectrum contains massless Goldstone bosons at infinite volume, this gives rise to 
finite size effects in finite volume quantities which vanish oc 1/L 2 . In ref. [30| this was found 
to hold for the scalar field expectation value, the Higgs mass and also for the fermion mass. 
In our model the Goldstone particles have a mass in an infinite volume due to the 0(4) 
symmetry breaking. On a finite volume this is expected to give rise to deviations from the 
linear 1/L 2 dependence only when the lattice extent L increases beyond the Goldstone cor- 
relation length, L > 0(1/ m^). Small deviations imply that the symmetry breaking effects 
are small. On the other hand, on the smaller lattices one expects additional non-leading 
finite-size effects. A pragmatic way to proceed is to apply the linear 1/L 2 extrapolation as 
long as no significant deviations are observed in the numerical data. 



In fig. H we plotted the observables vr, m a and tuf obtained at the three different points 
of region (I) as a function of 1/L 2 . Fig. |^ shows the 1/L 2 dependence of the same quantities 
for the points in the regions (II) and (III). The dashed lines are linear fits to all data points. 
It is seen that the linear 1/L 2 behavior holds satisfactory well within the statistical error 
bars in all three regions of the phase diagram. In region (I), where we expect the symmetry 
breaking effects to be even smaller than in the other regions, the agreement with a linear 
1/L 2 dependence is indeed best. In regions (II) and (III), only the mp values at the points 
which are most distant from the phase transition show a slight indication of the expected 
systematic deviation for large volumes. However, when comparing the results obtained in 
regions (II) and (III) to the ones of region (I), we conclude that this effect is rather small and 
we interpret this observation as further evidence that the effect of 0(4) symmetry breaking 
in our model is small. 



Also the couplings themselves show a 1/L 2 behavior as can be seen in fig. [10], where we 
have displayed as an example the ratio y/2X R = m a /vR (open symbols) as a function of 1/L 2 
for two points in region (II) and (III) in the phase diagram. We showed in the previous 
paragraph that both m a and Vr are proportional to 1/L 2 , vr ~ c\ + C2/L 2 , m a ~ di + d>2/L 2 , 
with relatively small coefficients C2 and d^. From an expansion in powers of 1/L 2 we ex- 
pect also the ratio m a /vR to obey an approximate 1/L 2 behavior. In addition we have also 



included in the graph the results for the corrected ratio y2X' R (filled symbols) which were 
obtained by means of eq. (4.11). Also here we expect, based on an expansion of the right- 
hand side of eq. ( 4.11|) in powers of 1/L 2 , an approximate linear 1/L 2 dependence which is 
also seen in the numerical results. Similar results hold also for the coupling y' R which we 
defined in eq. ( f4.10| ). In general we find that the linear 1/L 2 dependence is better fulfilled 
for the corrected ratios than for the uncorrected ones, which gives some evidence that the 
symmetry breaking effects in the corrected ratios are small. It is amusing that in all the 
cases the infinite volume extrapolations of the corrected and uncorrected ratios, shown by 
the dotted and dashed lines in fig. [TU], respectively, lead within the error bars to the same 



Figure 8: The quantities vr, m F and m a as a function of 1/L 2 for the three different 
coupling points of the region (I). The extent of the lattice in the time direction was kept fixed 
at T = 24. The dashed lines were obtained by a linear fit to all the data at different lattices 
and the symbols drawn at 1/L 2 = are the result of the infinite volume extrapolation. 



Figure 9: The same as fig. 8, but now for the regions (II) and (III). 



Figure 10: The ratios \/2\r = m a /vji (open symbols) and y2\' R (filled symbols) as a 
function ofl/L 2 for (n,y) = (0,4.0) (diamonds) and (—0.25,5.8) (squares). The extent in 
the time direction was kept fixed atT = 24. The dashed lines (for the open symbols) and the 
dotted lines (for the filled symbols) were obtained by linear fits to all data at different lattices 
and the symbols drawn at 1/L 2 = are the results of the infinite volume extrapolation. 





K 


y 




vr 


Vr 


m a /v R 


(I) 


0.29 
0.30 
0.31 


0.7 


0.259(8) 

0.62(2) 

0.81(2) 


0.08(2) 
0.21(2) 
0.271(2) 


1.00(9) 
0.87(3) 
0.866(8) 


2.8(5) 
2.9(2) 
2.96(8) 


(II) 


0.00 


3.8 
4.0 


0.41(4) 
0.601(7) 


0.1173(4) 
0.170(5) 


2.45(4) 
2.63(1) 


3.7(3) 
3.6(2) 


(III) 


-0.25 


5.6 
5.8 


0.46(4) 
0.56(4) 


0.114(2) 
0.150(2) 


2.56(5) 
2.7(1) 


4.1(3) 
3.8(3) 



Table 3: Infinite volume results for m a , vr, ur and m a /vR for the various points in the 
regions (I), (II) and (III). 



values, though the numerical values on the smaller volumes deviate by two or sometimes 
three standard deviations. 

In order to estimate errors we have performed extrapolations in various ways, e.g. we 
have used subsets of the lattices or we have either extrapolated the finite volume results of the 
couplings or computed their infinite volume estimates from the extrapolated infinite volume 
values of mjr, m a and Vr. The resulting variations are included in the errors quoted for the 
infinite volume results. In table || we give the infinite volume results for the renormalized 
couplings. An estimation of the symmetry breaking corrections as discussed in sect. 4 shows 
that these effects are small, less than 6% in a parameter range with mp < 0.5 and m a < 0.7. 

5.2 Renormalized couplings 

A comprehensive summary of our infinite volume results for the renormalized couplings is 
given in fig. [l| where we have plotted the ratio m a jvR, = ^/2\r as a function of mp/vR = ur 
for the various points of the regions (I), (II) and (III). The symbols were chosen such that 
they match with those in fig. § so that the reader can easily find out at which bare couplings 
in the phase diagram the various points have been obtained. All points have roughly the 
same value of the cut-off in units of the scalar field expectation value, clvr ~ 0.15-0.25. 
The unsystematic scattering (as a function of the cut-off) of the points in the region around 
Ur = 2.5, \/2Ar = 4 is still within one standard deviation and is presumably due to uncer- 
tainties in the infinite volume extrapolation. The figure shows that the data points obtained 
in the regions (II) and (III) fall almost on top of each other. The renormalized couplings 
appear to saturate when k \ and not to increase further when k is lowered beyond zero, 
always keeping av r roughly fixed. 

We furthermore compare in fig. [I] the numerical results for a/2A^ and jjr with the tree 
level unitarity bounds marked by the horizontal and vertical arrows below and beside the 
axis. The tree level unitarity bound for a/2Ar was taken from ref . || . The bound on vr for 
a fermion-Higgs model with No doublets is given by 

y R <^Aix/N D , (5.1) 

i. e. yR^j2.h for Np, = 2. The points obtained in the regions (II) and (III) are very close 
to the tree level unitarity bounds which indicates that the couplings are not very strong. 



This is also consistent with our observation that the numerical results for the scalar propaga- 
tors could well be described in the full momentum range by the one fermion loop expressions. 



For comparison we have displayed in fig. [I] also the yR dependence of \/2\ R using the 
assumption that the one-loop /3-function is a good approximation to the full /3-function. The 
one-loop /3-functions for the Yukawa and quartic self-coupling in a fermion-Higgs model with 
Ne> doublets are given by the expressions: 

CM = (5-2) 



dy(t) 


= ND f 


dt 


4tt 2 ^ 


dX(t) 


- 3 A 2 


dt 


" 2tt» A 



>'a07.A) = ^ L = ^^ + ^\y 2 -^t , (5-3) 



7T^ 7T^ 



where A and y are the running couplings and t = ln(/x/Ai). The quantity fi is an energy 
scale and A x is the cut-off of the one-loop /3-function model. In order to find the relations 
for the renormalized couplings we have integrated these /3-functions numerically from the 
cut-off scale \i = Ai down to the physical scale \i = Vr, with the identifications y(0) = y , 
Vita) = Vr, where tn = ln(t>^/Ai). Since the simulations were carried out at A = oo, we have 
chosen Ao = 100 ~ oo and varied the starting value for yo in the interval < y < 50. The 
solid and dotted curves in fig. [L] correspond to two different values of the ratio Vr/ A\. For 
the solid curve this ratio was adjusted such that the agreement with the numerical results is 
best. It is remarkable that the shape of the curve is in reasonable agreement with our data. 
The ratio vr/ Ai of the dotted curve is by a factor two smaller than that of the solid curve. 

Let us now turn to the question of the determination of upper bounds on m a and mp- 
Reading off an upper bound at m a = 0.7/a, we find m cr /vji'^3 in region (I) of the phase 
diagram and m a /vR^A in regions (II) and (III). The fermion mass increases along the phase 
transition line (keeping am a constant) until y ~ 4 where K c (y) becomes negative. From then 
on it appears to remain roughly constant. Also the field expectation value does not change 
much for y^A and we find wir/vr = 0.9, 2.6 and 2.6 in regions (I), (II) and (III) respectively. 
The observed saturation of the renormalized Yukawa coupling indicates that we can read off 
the upper bound for y R at y ~ 4. 

It is interesting to compare our results for Xr with those obtained in the mirror fermion 
model of ref. fl9fl , which also describes two doublets of light fermions. Since the mirror 



fermion model has a quite different bare action with a different bare parameter space, a 
comparison of the yR dependence of Xr at the same values of civr would give some idea to 



what extent these results are model dependent. Such a comparison is shown in fig. |TT| where 
we have plotted the ratio m a /vR against ttir/vr for various models. The three diamonds 
represent some infinite volume results in the reduced staggered fermion model which have 
roughly the same cut-off, avR = 0.15-0.21. In contrast to our results, the results in the 
mirror fermion model (full circles) |2S| were not extrapolated to infinite volume. The cut-off 
values on the 6 3 12 lattice, avR = 0.31-0.39, are consistent with our results on the 6 3 24 lattice 



see table 2). It can be seen from fig. [11] that the points obtained from the mirror model are 



within error bars on the same curve as our results. 



We included in fig. |TT] also some infinite volume results which were obtained in the pure 
0(4) model (N D = 0, square, avR = 0.20(2)) |35j and in the naive fermion model (N D = 32, 
open circle, qvr = 0.20(3)) PU|. It is seen that in the naive fermion model, where uif/vr is 



relatively small, the fermion effect on Xr is negligible, whereas the fermions give a moderate 



Figure 11: Comparison of the results for the ratio m a /vR as a function of mp/vR obtained 
in various models. The open symbols represent infinite volume results for the 0(4) model 
(Nd = ; square) fiSh]] , the naive fermion model (Njj = 32, open circle) and the 

reduced staggered fermion model (Njj = 2, diamond), whereas the full circles show finite 
volume (6 3 12) results obtained in the mirror fermion model (Np = 2) fi2b]] . 



increase of Xr in the model with two doublets. In both models the maximal values for 
appear to be close to the tree level unitarity bound given in eq. (|5.1|) . 



An interesting speculation has been that, even though the Yukawa and quartic self- 
couplings are trivial, they might still be sufficiently large as to give rise to interesting non- 
perturbative effects, like formation of a p bound state We find, however, that the 

largest values of Xr and dr we can obtain are very close to their tree level unitarity bounds, 
which suggests that they are still on the edge of the perturbative regime. This indication is 
supported by the observation here and also in ref. |30| for the 32 doublet model, that the 



numerical results for the scalar propagator can be very accurately described using one-loop 
renormalized perturbation theory. It is therefore most likely that renormalized perturbation 
theory gives a complete and accurate description of the physics of Yukawa models. 



6 Summary and discussion 

In this paper we have used 'reduced' or 'real' staggered fermions in a high statistics investi- 
gation of a fermion-Higgs model. The two staggered flavors are coupled to the 0(4) Higgs 
field, leading to a model with a single isospin doublet in the scaling region. In a simulation 
with the Hybrid Monte Carlo algorithm this number of doublets has to be doubled. Since 
the species doubler degrees of freedom are used as physical flavor-spin components, there is 
no redundancy in the fermion field in this approach and the model can be very efficiently 
simulated on large lattices. 

The Yukawa coupling to the staggered flavors breaks the 0(4) symmetry and it requires 
two scalar field counterterms to restore this symmetry in the scaling region. Even without 
these counterterms, this 0(4) symmetry breaking has only a small effect on the values of 
the renormalized couplings, at least in the physically relevant scaling region of the phase 
diagram where we have done our simulations. This follows from our findings that a) the 
deviations from the 1/L 2 finite size effects are very small, b) the Goldstone mass is much 
smaller than the Higgs mass and c) that the one fermion loop estimate of the symmetry 
breaking correction to the quartic coupling is found to be small. The one fermion loop esti- 
mate for the symmetry breaking effects is presumably reliable, because is gives a reasonable 
prediction for the Goldstone mass. 

The renormalized Yukawa coupling and the renormalized quartic self-coupling Xr were 
determined from the tree level relations yn = ttif/vr and Xr = m^/2v R , where vr = v / 
is the renormalized field expectation value. On a finite lattice the masses mp and m a and 
the Goldstone wave-function renormalization constant Z v were computed from the scalar 
and fermion propagators in momentum space. The momentum dependence of the scalar 
propagators differs from that of a free propagator and in particular on small lattices it is 
important to include the one fermion loop self-energy to get a reliable estimate for m a and 
Z^. In order to investigate the finite size dependence of our results and extrapolate them 
to the infinite volume, we have simulated on a sequence of lattices ranging in size from 6 3 24 
to 16 3 24. The finite volume results for m a , vr and tuf show a 1/L 2 dependence, which was 
used to extrapolate these observables to the infinite volume. 

A comprehensive summary of our infinite volume results is shown in fig. [I], where we have 
plotted the ratio m a /vR = \/2Xr as a function of tjif/vr = yR. The results were obtained 



in three different regions of the phase diagram (c.f. fig. ||]) and always at A = oo. The ratio 
m a/vR increases a little bit for increasing jjr and becomes at large i/r slightly larger than the 
tree level unitarity bound. The maximal value of the Yukawa coupling is roughly equal to 
its tree level unitarity bound. This indicates that the couplings are not very strong. Further 
evidence for the perturbative nature of the renormalized couplings comes from the fact that 
the numerical results for the scalar field propagators can be nicely described by taking into 
account the one fermion loop contribution to the self-energy. We find even for the points in 
region (III) an almost perfect agreement over the full range of momenta. Furthermore fig. [I] 
shows that the numerical values for m a /vR and tuf/vr are in qualitative agreement with a 
simple model based on extending the one- loop /3-function to infinite couplings (solid curve). 
In conclusion we can say that all our results are consistent with the triviality scenario and 
there is no evidence for a strong coupling sector. 



Our results are obtained in a particular regularization of the SU(2) fermion-Higgs model: 
using the lattice regularization with reduced staggered fermions and two isospin doublets in 
the scaling region. One would like to see to what extent the results are model dependent. 
Therefore we have compared our results with those of ref. |25| which uses mirror fermions 



with two doublets in the scaling region. Looking at the i/r dependence of Xr at roughly the 
same value of the cut-off in all models, we find consistent results within error bars, c.f. fig. |TT|. 

The main result of our investigation therefore is that the renormalized quartic and Yukawa 
couplings cannot be strong, unless the cut-off is unacceptably low. In fact we find that one- 
loop renormalized perturbation theory is applicable even for the maximally strong couplings. 
More quantitatively we find for the upper bounds on the masses of the Higgs particle and 
the heavy fermion in our model, m a ^AvR and mF^2.6vR for cut-off values A^2.5 m aF . 
From experience in the 0(4) model with modifications in the regularization (by including 
dimension six operators [JTOj or using improved actions Jll|]) we expect that these numbers for 
the upper bounds may be stretched by perhaps 20-30%. Larger changes have been recently 



found with the Pauli-Villars regularization scheme in the large N limit [13 
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